Two-dimensional array of magnetic particles: The role of an interaction cutoff 
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Based on theoretical results and simulations, in two-dimensional arrangements of a dense dipolar 
particle system, there are two relevant local dipole arrangements: (1) a ferromagnetic state with 
dipoles organized in a triangular lattice, and (2) an anti-ferromagnetic state with dipoles organized in 
a square lattice. In order to accelerate simulation algorithms we search for the possibility of cutting 
off the interaction potential. Simulations on a dipolar two-line system lead to the observation that 
the ferromagnetic state is much more sensitive to the interaction cutoff R than the corresponding 
anti-ferromagnetic state. For R > 8 (measured in particle diameters) there is no substantial change 
in the energetical balance of the ferromagnetic and anti-ferromagnetic state and the ferromagnetic 
state slightly dominates over the anti-ferromagnetic state, while the situation is changed rapidly 
for lower interaction cutoff values, leading to the disappearance of the ferromagnetic ground state. 
We studied the effect of bending ferromagnetic and anti-ferromagnetic two-line systems and we 
observed that the cutoff has a major impact on the energetical balance of the ferromagnetic and 
anti-ferromagnetic state for R < 4. Based on our results we argue that R « 5 is a reasonable 
choice for dipole-dipole interaction cutoff in two-dimensional dipolar hard sphere systems, if one is 
interested in local ordering. 
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I. INTRODUCTION 

Long range interaction represents a major challenge for 
computer simulations. The size of the tractable systems 
(N particles) is limited through the fact that the order of 
N 2 calculations are to be carried out at each step, though 
for many purposes large systems need to be studied. Pe- 
riodic boundary conditions, which are often helpful, can 
be implemented only by using sophisticated summation 
algorithms (if possible due to screening) . 

In principle the above problems occur in so called short 
range interaction models as well, like in the most exten- 
sively studied Lennard-Jones system. However, for these 
systems a cutoff is usually introduced making the original 
short range model explicitly finite range. It is generally 
accepted that the error introduced by the cutoff is negli- 
gible provided the cutoff distance is large enough 

Frequently the long range interaction potential falls off 
like r -1 where r is the distance between the particles. 
This should be compared to the Lennard-Jones system 
where the potential decreases like r -6 , where it is as- 
sumed that the attractive part of that potential is due to 
induced dipole-dipole interaction. 

In this study we focus on the question of cutting off a 
potential which is in between the two above cases. We 
consider a two-dimensional ensemble of magnetic parti- 
cles interacting with a r~ 3 potential, which, however, has 
an orientation dependence as well. It is crucial from the 



point of view of efficient programing to know if a reason- 
able cutoff can be introduced in this system. We inves- 
tigate this problem by comparing the stability of static 
configurations. 

II. THE LUTTINGER-TISZA METHOD 

The Hamiltonian of a system of spherical dipoles is 

H=lj2sJj ijSj , (1) 

where i and j are dipole indexes, s denotes the dipole 
momentum vector, s T denotes the transpose of s, and 

1 / Yij O Yij \ 

where I denotes the identity matrix, and r.^ denotes the 
relative position vector of two dipoles. The 1/2 factor in 
Eq. avoids double counting of dipole pairs. 

We can study the crystalline state of a dipole system 
using the Luttinger-Tisza method [2j based on the idea 
that in case of crystals it is a natural assumption that the 
ground state exhibits some discrete translational sym- 
metry. If r(i) denotes the points generated from i with 
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discrete translations belonging to the T symmetry group 
the mentioned symmetry corresponds to Sj = s,/ for all 
i' S r(i). According to this the system can be broken 
into identical cells and the summation in Eq. can be 
limited to summation over one single cell. According to 
this the energy per dipole can be expressed as 
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where n is the number of dipoles per cell, and Ay are 
symmetric matrices defined by 
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The expression of the energy per dipole in Eq. can 
be simplified by considering the s = (si)" =1 hyper- vector 
and A = (Ajj)™,' =1 hyper-matrix. By construction A is 
symmetric. Using these, E can be written as 
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Solving the eigenvalue problem of the nd dimensional 
symmetric matrix A, where d is the dimension of the 
dipoles, we find the A^ eigenvalues and orthogonal 
eigenvector system with normalization ||xfe|| = y/n. Us- 
ing these we have 
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where bk denotes the components of s in the x^ orthog- 
onal eigenvector system. 

If the dipoles have identical scalar strength fi, i.e. 
||sj|| = fi, then bk must satisfy for all i = l...n the 
condition 
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where x % k are the components of the x^ hyper-vector 
which belong to dipole index i. Adding the square of 
the above equations and taking into consideration that 
{kk}k=i form an orthogonal system, we conclude that bk 
must satisfy the condition 
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In the framework of the Luttinger-Tisza method these 
two conditions are known as the strong (Eq. Q) and the 



weak (Eq. (JSJ) conditions. From the weak condition and 
Eq. © it can be derived that the energy per dipole in the 



ground state is E. n 



1/2 A, 



2 , where X m in denotes 



the smallest eigenvalue of A. If there is one single eigen- 
value equal to A TO j n the ground state dipole arrangement 
is given by the corresponding eigenvector. If there are 
more eigenvalues equal to \ m im the ground state dipole 
arrangement is given by the linear combinations of the 
corresponding eigenvectors which satisfy the strong con- 
dition. 



III. TWO-DIMENSIONAL ARRAY OF 
MAGNETIC PARTICLES 

The above method was applied to a system of two-di- 
mensional dipole moments with identical scalar strength 
located at the sites of an infinite rhombic lattice with 
an arbitrary rhombicity angle by Brankov and Danchev 
[3j . They considered that the ground state of this system 
has a translational symmetry corresponding to discrete 
translations along the lattice lines with multiples of 2a, 
where a is the lattice constant. They found that the 
ground state depends on the rhombicity angle. 

We repeated their calculations with the consideration 
that the dipoles are carried by identical hard spheri- 
cal particles of diameter equal to the lattice constant, 
and according to the geometrical constraint introduced 
by this consideration we limited the rhombicity angle to 
60° < a < 90°. In accordance with their results we 
found that the system has a ferromagnetic ground state 
for 60° < a < 79.38°, and an anti- ferromagnetic ground 
state for 79.38° < a < 90°. We also found that the 
ground state for a = 60° is a continuously degenerate 
ferromagnetic state, and for a — 90° is a continuously de- 
generate microvortex state including a four-fold degener- 
ate anti-ferromagnetic state, where the microvortex state 
is defined as two anti-ferromagnetic sublattices making 
an arbitrary angle with each other. We also identified 
the states with the second lowest energy per dipole. The 
results for a = 60° and a = 90° are summarized in Tab. [I] 

We repeated the calculations taking into consideration 
the interaction of only two neighboring lines on the rhom- 
bic lattice. This corresponds to the interaction of two 
lines of dipolar hard spheres shifted according to the a 
rhombicity angle. The Luttinger-Tisza method can be 
applied in a straightforward way also in this case. We 
observed that the ground state depends on the rhom- 
bicity angle a similar to the previous case. We found 
that the system has a ferromagnetic ground state for 
60° < a < 75.67°, and an anti- ferromagnetic ground 
state for 75.67° < a < 90°. The ground state for 
a = 60° is a two-fold degenerate ferromagnetic state, and 
for a — 90° is a two-fold degenerate anti-ferromagnetic 
state. We also identified the states with the second low- 
est energy per dipole, and we summarized the results for 
a = 60° and a = 90° in Tab. HJ 

It is not surprising that taking into consideration only 
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a = 60° 


E = -2.758 


continuously degenerate 
ferromagnetic state 


E = -2.047 


six-fold degenerate 
anti-ferromagnetic state 


a = 90° 


E = -2.549 


continuously degenerate 
microvortex state including 
a four-fold degenerate 
anti-ferromagnetic state 


E = -2.258 


continuously degenerate 
ferromagnetic state 



TABLE I: Dipole arrangements corresponding to the lowest 
two energy values per dipole in a two-dimensional system of 
dipole moments with identical scalar strength located at the 
sites of an infinite rhombic lattice with rhombicity angle a = 
60° and a = 90°. The energy is measured in /i 2 /a 3 units, 
where fi is the scalar strength of the dipole moments and a is 
the lattice constant. 



a = 60° 


E = -2.582 


two-fold degenerate 
ferromagnetic state 


E = -2.226 


two-fold degenerate 
anti-ferromagnetic state 


a = 90° 


E = -2.477 


two-fold degenerate 
anti-ferromagnetic state 


E = -2.331 


two-fold degenerate 
ferromagnetic state 



TABLE II: Dipole arrangements corresponding to the lowest 
two energy values per dipole in a two-dimensional system of 
dipole moments with identical scalar strength located on two 
neighboring lines of an infinite rhombic lattice with rhombic- 
ity angle a = 60° and a = 90°. The energy is measured 
in fJ, 2 /a 3 units, where fj, is the scalar strength of the dipole 
moments and a is the lattice constant. 



two lines of the rhombic lattice reduces significantly the 
original symmetry of the system. This can be seen com- 
paring the results in Tab.[I]and Tab.lTTl It is important to 
note that the two-line system has no continuously degen- 
erate ground state, and thus the ground state is always 
defined by the Luttinger-Tisza basic arrangement with 
the lowest eigenvalue. In particular there is a special 
ferromagnetic and an anti-ferromagnetic state, which - 
depending on the a rhombicity angle - define the ground 
state. 

It can be also seen comparing the results in Tab. [I] and 
Tab. El that the interaction of two neighboring lines al- 
most saturates the long-range dipole-dipole interaction, 
furthermore it is widely known that dipolar spheres due 
to dipole-dipole interactions tend to aggregate into chain 
like structures (see for example 0, ]f| and references 
therein) in which the energies of intrachain interactions 
are much greater than those of interchain interactions 



[6j. These confirm that studying a two- line system gives 
valuable results related to properties of dipole-dipole in- 
teraction in general. 



IV. THE ROLE OF AN INTERACTION 
CUTOFF 

Brankov and Danchev observed that the ground 
state of a system of dipoles on an infinite rhombic lattice 
is sensitive to the dipole-dipole interaction range. An R 
interaction cutoff distance can be introduced in a natural 
way with a slight modification of Eq. Q as 
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where denotes the distance of two dipoles, and the 
Tif < R constraint represents the fact that the summa- 
tion must only contain terms corresponding to pairs of 
spherical dipoles closer to each other than the R cutoff 
distance. We measure the interaction cutoff value in a 
units, where a is the lattice constant equal to the parti- 
cle diameter. 

It can be seen from Eq. @ that the strength of the 
dipole-dipole interaction decays with 1/rfji, and thus the 
above expression for large R can be arbitrarily close to 
the long-range limit in Eq. (0}. This means that the 
numerical evaluation of Luttinger-Tisza states in general 
can be based on Eq. © if R is big enough. 

Our numerical results for the two-line system at R 
equal to 10 6 (lines) and 8 (points) are shown on Fig. ^ 
The results corresponding to R = 10 6 are close to the 
long-range interaction limit within the numerical errors 
of 64 bit floating point arithmetic. The upper panel of 
Fig. ^ shows the lowest energy per dipole of the ferro- 
magnetic and the anti-ferromagnetic state as function of 
the rhombicity angle. The lower panel of Fig. ^ shows 
the angle of the dipoles with respect to the direction of 
the longest linear dimension of the system. 

Without any calculation one might expect that in the 
ground state of the two-line system the dipoles are ori- 
ented parallel to the lines in both ferromagnetic and anti- 
ferromagnetic state. It is a surprising result of our cal- 
culations that this is true only for a = 60° and a = 90°. 
For any other a the dipoles form a small angle with the 
lines. However these angles are less than 2.5° they cannot 
be neglected. 

Below a certain a, as can be observed on Fig. ^ the 
ground state of the system corresponds to the ferromag- 
netic order, and above it to the anti-ferromagnetic order. 
The angle at which the transition from a ferromagnetic to 
an anti-ferromagnetic ground state happens, is shifted by 
only 3% due to the cut off. That the anti-ferromagnetic 
state remains almost unchanged is a consequence of the 
strong coupling of neighboring dipoles of opposite orien- 
tation, which makes the interaction cutoff irrelevant. 
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FIG. 1: Numerical results for the two-line system at R equal 
to 10 6 (lines) and 8 (points). The upper panel shows the 
lowest energy per dipole of the ferromagnetic and the anti- 
ferromagnetic state as function of the rhombicity angle. The 
energy is measured in units fi 2 /a 3 . The lower panel shows the 
angle in degrees which the dipoles form with the direction of 
the longest linear dimension of the system. 

As expected (see Fig. ^) the ferromagnetic state at 
a = 60° and the anti-ferromagnetic state at a = 90° 
are stable configurations of the two- line system of dipolar 
hard spheres independent of R, however as R decreases 
the crossover point gets slightly shifted. 

At low interaction cutoff distances (R < 8) the discrete 
nature of the system becomes more and more relevant 
and both the ferromagnetic and anti-ferromagnetic en- 
ergy per dipole begin to exhibit sudden jumps in function 
of the rhombicity angle (see Fig. 0). As the interaction 
cutoff decreases the energy jumps become more and more 
relevant. This behavior can introduce numerical instabil- 
ities in simulations using a badly chosen cutoff distance. 
It must be noted however that one should not overesti- 
mate this effect as the energy jumps are relatively small. 
It is an interesting observation that for some interaction 
cutoff values (e.g. at R « 4) the energy per dipole shows 
significantly lower anomalies. 

The ferromagnetic line is shifted upward as R decreases 




Rhombicity angle (a) 




Rhombicity angle (a) 



FIG. 2: Numerical results for the two-line system at R equal 
to 7, 6, 5, 4, and 3. The upper panel shows the lowest energy 
per dipole of the ferromagnetic and the anti-ferromagnetic 
state as function of the rhombicity angle. In ferromagnetic 
states the lines are shifted upward as R decreases. The lower 
panel shows the angle in degrees which the dipoles form with 
the direction of the longest linear dimension of the system. 



(see Fig. ^ and Fig. and according to this the anti- 
ferromagnetic state becomes more and more dominant. 
For large R the ferromagnetic state at a = 60° has lower 
energy per dipole than the anti-ferromagnetic state at 
a = 90°. Our numerical results show that at R « 4 
the situation is reversed, and at R w 2 the ferromagnetic 
ground state disappears. Brankov and Danchev Q found 
that in case of an infinite rhombic lattice with rhombicity 
angle a = 60° the ferromagnetic ground state disappears 
at R w 3. 



V. FINITE SIZE CORRECTIONS 

We investigated the finite size corrections of the en- 
ergy per dipole of the two-line system in ferromagnetic 
and anti-ferromagnetic states. In these states the infinite 
system can be decomposed into identical finite segments. 
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If N denotes the number of dipoles per line in a finite 
segment, the energy per dipolc of the infinite system can 
be written as 
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where a(N) denotes the dipoles belonging to one seg- 
ment, and <t(N) c denotes the complementer of <j(N). 
The first part in the above expression can be recognized 
as the energy £(N) per dipole of a finite segment con- 
taining N dipoles per line. 

We define the following quantity of energy dimension 



d£{N) = N £(N) - E 
It can be seen that 
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where in case of an interaction cutoff R one may add the 



condition 



7*jj < R. As is proportional to 1/rfj (see 
Eq. J2J), one may expect that for large system size d£(N) 
is independent of N, and thus the limit 



dE= lim d£(N), 



(13) 



exists and is finite. Our numerical investigations con- 
firmed this expectation. The convergence of d£(N) is 
of order 1/N in the ferromagnetic case, and is of order 
1/N 3 in the anti-ferromagnetic case. This proves that 
the above quantity is well defined. We refer to the above 
quantity as the finite size coefficient. 

Numerical results showing the dependence of the en- 
ergy per dipole £(N) of the two-line system on the sys- 
tem size in the long-range limit are presented in Fig. [3] 
The upper panel shows results for N equal to 10 5 (lines) 
and 100 (points), and the lower panel shows results for 
N equal to 10, 8, 7, 6, and 5. Both the ferromagnetic 
and anti-ferromagnetic lines are moved upward as N de- 
creases. £(N) at N = 10 5 is close to the energy per 
dipole E of the infinite two-line system within the nu- 
merical errors of 64 bit floating point arithmetic. 

Based on the definition of the finite size coefficient for 
large N the energy per dipole of a finite system can be 
approximated as 



£{N) « E + dE/N. 



(14) 
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FIG. 3: Numerical results showing the dependence of the 
energy per dipole of the two-line system on the system size 
in the long-range limit. The upper panel shows results for 
N (number of particles per line) equal to 10 5 (lines) and 100 
(points). The lower panel shows results for equal to 10, 8, 
7, 6, and 5. Both the ferromagnetic and anti-ferromagnetic 
lines are moved upward as JV decreases. (Note the different 
scales on the vertical axes.) 



Our numerical investigations show that this approxi- 
mation is reasonable even for N m 10. The finite size 
coefficient of the ferromagnetic state is approximately 
two times bigger than the finite size coefficient of the 
anti-ferromagnetic state, and thus the ferromagnetic line 
moves upward approximately two times faster than the 
anti- ferromagnetic line (see Fig. It can be observed 
that for large N the ferromagnetic state at a = 60° has 
lower energy per dipole than the anti-ferromagnetic state 
at a = 90°. Our numerical results show that at N = 20 
the situation is reversed, and at N = 5 the ferromagnetic 
ground state disappears. 
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FIG. 4: Numerical results showing the dependence of 
the finite size coefficient for both ferromagnetic and anti- 
ferromagnetic state at R — I0 6 (lines), for ferromagnetic 
state at R equal to 1000 and 500 (points), and for anti- 
ferromagnetic state at R equal to 100 and 50 (points). The 
upper panel shows results for the ferromagnetic state and 
the lower panel shows results for the anti-ferromagnetic state. 
The finite size coefficient is measured in units (i 2 /a 3 . In the 
ferromagnetic case the lines are lowered as R decreases. (Note 
the different scales on the vertical axes.) 



VI. DEPENDENCE OF THE FINITE SIZE 
COEFFICIENT ON INTERACTION CUTOFF 

We investigated the dependence of the finite size co- 
efficient on the interaction cutoff distance R. Fig. 0] 
shows the numerical results for ferromagnetic state (up- 
per panel) at R equal to 1000 and 500 (points), for anti- 
ferromagnetic state (lower panel) at R equal to 100 and 
50 (points), and for both ferromagnetic (upper panel) 
and anti- ferromagnetic state (lower panel) at R = 10 6 
(lines). The finite size coefficient is measured in units 
fi 2 /a 3 . We calculated its value by evaluating the expres- 
sion in Eq. JEJ at N = 10 5 . The results for R = 10 6 are 
close to the long-range limit within the errors of 64 bit 
floating point arithmetic. 

As R is lowered in the anti-ferromagnetic case the fi- 



nite size coefficient remains almost unchanged even for 
R w 50, while in the ferromagnetic case it decreases sig- 
nificantly already at R w 1000. This shows again that 
the ferromagnetic state is much more sensitive to the in- 
teraction cutoff than the anti-ferromagnetic state. 

At lower interaction cutoff distances (at R < 50) the 
discrete nature of the system manifests itself in sudden 
jumps in the finite size coefficient (see Fig- EJ . 

The upper panel of Fig. [3] shows the finite size coeffi- 
cient of the ferromagnetic state as function of the rhom- 
bicity angle at R equal to 4, 3.75, 3.5, 3.25, and 3. The 
lines are lowered as R decreases. The lower panel shows 
the finite size coefficient of the anti-ferromagnetic state 
at R equal to 40, 16, 8, 4, and 3. The lines are shifted 
upward as R decreases. 

The jumps in the finite size coefficient become bigger as 
the interaction cutoff decreases (see Fig.|SJ. These jumps 
are not relevant at large N, but can introduce energy 
jumps at lower dipole numbers, and thus can introduce 
local numerical instabilities in simulations, but this effect 
should not be overestimated as the introduced energy 
jumps are relatively small. 



VII. BENDING TWO LINES OF MAGNETIC 
PARTICLES 



The finite size behavior presented before gives a good 
description of finite dipole systems at large N, but it is 
not too helpful at lower N. For a better understand- 
ing of the system, we studied numerically finite sys- 
tems (at small N) investigating the effect of bending two 
lines of dipolar hard spheres in ferromagnetic and anti- 
ferromagnetic states (see Fig. El (a) and (b)). In unbent 
case these correspond to the previously studied ferromag- 
netic state at a — 60° and anti-ferromagnetic state at 
a = 90°. We introduce the 7 bending parameter and we 
define the bent system as composed of particles placed on 
an arc of angle 2Nj with dipole vectors tangential to the 
arc (see definition of 7 on Fig. 0(a) and (b)). This defini- 
tion involves a so called 'bending limit' as the arc's angle 
is limited to 2n, and thus 7 must satisfy the 7 < n/N 
condition. 

Our numerical results show that for bending either a 
ferromagnetic or anti-ferromagnetic two- line system some 
physical effort is needed. We observed that the two-line 
system in ferromagnetic state can be bent easier than 
in the corresponding anti-ferromagnetic state. This is a 
consequence of the strong coupling of neighboring dipolcs 
oriented anti-parallel. We also observed that as the anti- 
ferromagnetic state is bent it becomes less and less sta- 
ble. Fig. shows numerical results related to bending 
a two-line system at different interaction cutoff values. 
As function of the bending parameter 7 and system size 
N we compared the energy per dipole of the ferromag- 
netic and anti-ferromagnetic states and we identified the 
points (7, N) at which these two states are energetically 
equivalent. We repeated this procedure at different R in- 
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FIG. 5: Numerical results showing the dependence of the fi- 
nite size coefficient at lower interaction cutoff distances. The 
upper panel shows the finite size coefficient of the ferromag- 
netic state as function of the rhombicity angle at R equal 
to 4, 3.75, 3.5, 3.25, and 3. The lines are lowered as R de- 
creases. The lower panel shows the the finite size coefficient 
of the anti-ferromagnetic state at R equal to 40, 16, 8, 4, and 
3. The lines are shifted upward as R decreases. (Note the 
different scales on the vertical axes.) 



FIG. 6: Numerical results related to bending a two-line sys- 
tem at different interaction cutoff values. The upper panel 
shows a finite system of two lines of dipolar hard spheres in 
ferromagnetic and anti-ferromagnetic states. The lower panel 
shows (7, N) state diagrams (see text for description) for R 
ranging from 2 to 00. The lines are moved upward and lower 
as R decreases. 



VIII. CONCLUSIONS 



teraction cutoff values. The lower panel of Fig. shows 
corresponding (7, N) state diagrams. 

In the long-range limit for small system size and low 
bending parameter the anti-ferromagnetic state has lower 
energy per dipole. This is in accordance with our pre- 
vious results, and remains valid up to R « 5. It is a 
surprising result that this behavior changes rapidly for 
interaction cutoff values between 4 and 5. For R < 4 the 
anti-ferromagnetic state remains more stable at large N 
values even for large bending parameters. This means 
that at this point the general characteristics of an arbi- 
trary dipole system is substantially changed. Based on 
Fig. [5] and on our previous results we argue that R ~ 5 
is a reasonable choice for dipole-dipole interaction cutoff 
for two dimensional systems of dipolar hard spheres, if 
one is interested in local ordering. 



Based on the fact that dipolar spheres due to dipole- 
dipole interactions tend to aggregate into chain like struc- 
tures in which the ratio of interchain-to-intrachain inter- 
actions is small, and that moreover the interaction of par- 
allel chains of dipolar hard spheres almost saturates the 
dipole-dipole interaction in two-dimensional dense sys- 
tems, we argue that the study of a dipolar two-line system 
gives valuable results for general dipolar particle systems. 

Theoretical results and simulations show two relevant 
dipole arrangements: (1) a ferromagnetic state with 
dipoles organized in a triangular lattice, and (2) an anti- 
ferromagnetic state with dipoles organized in a square 
lattice. Numerical results on a dipolar two-line system 
show that the ferromagnetic state is much more sen- 
sitive to the interaction cutoff than the corresponding 
anti-ferromagnetic state. This can be explained by the 
efficient coupling of dipoles oriented anti-parallel. For 
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R > 8 there is no substantial change in the energetical 
balance of the ferromagnetic and anti- ferromagnetic state 
and the ferromagnetic state slightly dominates over the 
anti-ferromagnetic state, while the situation is changed 
rapidly for lower interaction cutoff values, leading to the 
disappearance of the ferromagnetic ground state. Our 
numerical results show that the ferromagnetic ground 
state disappears at R w 2. Brankov and Danchev |3J 
found that in case of an infinite triangular lattice the 
ferromagnetic ground state disappears at R 3. 

For characterizing the finite size behavior of the two- 
line system we introduced a finite size coefficient, and 
we observed that it is sensitive to the interaction cut- 
off for both ferromagnetic and anti- ferromagnetic states. 
We also observed that at low interaction cutoff values the 
discrete nature of the system leads to small energetical 
anomalies. These anomalies increase as the interaction 
cutoff is lowered and can introduce instabilities in nu- 
merical simulations. We argue however that this effect 
becomes relevant only at first or second neighbor inter- 
action and it can be neglected at higher interaction cutoff 
values. 

Finally we studied the effects of bending ferromagnetic 
and anti-ferromagnetic two-line systems. We character- 
ized the bending of a two-line system with the parameter 
7, while N is the number of dipoles per line. We cre- 
ated (7, N) state diagrams separating energetically favor- 
able ferromagnetic and anti- ferromagnetic states. We ob- 
served that there is a substantial change of these state di- 
agrams for R < 4, and - in accordance with our previous 
results - we argue that R w 5 is a reasonable choice for 
dipole-dipole interaction cutoff in two-dimensional dipo- 
lar hard sphere systems, if one is interested in local or- 
dering. 

It is a surprising result that the reasonable interaction 
cutoff is independent of the strength of the dipole-dipole 
interaction and the particle size. This is a consequence of 
the fact that there are two relevant dipole arrangements 
(a ferromagnetic and an anti- ferromagnetic), and their 
energetical balance can be reduced to geometrical factors. 
If there are any other interactions in the system (e.g. 
friction), this study must be revisited and it may turn 
out that the reasonable interaction cutoff is dependent 



on the interaction strength and particle size. We envision 
however that in some cases (e.g. in case of friction) the 
presence of another short-range interaction keeps or even 
lowers the value of the reasonable interaction cutoff found 
above. 

In this paper we focused on the local dipole order- 
ing. In the ferromagnetic case however domain structures 
become important, which can reduce external magnetic 
stray fields. These global structures should depend on the 
long range part of the interaction. For magnetic granular 
systems the formation of such domains may be hindered 
e.g. by friction, though, as it requires the reorientation 
of particles. 

We did not address the response to an external mag- 
netic field. The reason is that long range correlations 
and hence the response functions will be more strongly af- 
fected by a dipolar interaction cutoff than the local struc- 
tures and energy densities considered in this paper. In 
principle an Ewald summation method 0, Q would allow 
to explore the response properties in the thermodynamic 
limit in terms of large but finite systems with periodic 
boundary conditions. However, here again friction may 
be an important factor to be taken into account: An 
external magnetic field trying to orient the magnetic mo- 
ments would exert a stress on the particle arrangement, 
if particle rotations would be hindered by friction. Then 
the magnetic response of the system would crucially de- 
pend on the relative strength of the magnetic anisotropy 
of the particles, i.e. the coupling between particle and 
magnetic moment orientations, and friction forces be- 
tween the particles. 
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